Magnetization and specific heat of the dimer system CuTe 2 Os 
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We report on magnetization and specific heat measurements on single-crystalline CuTe2 0s. The 
experimental data are directly compared to theoretical results for two different spin structures, 
namely an alternating spin-chain and a two-dimensional (2D) coupled dimer model, obtained by 
Das et al. [Phys. Rev. B 77, 224437 (2008)]. While the analysis of the specific heat does not allow 
to distinguish between the two models, the magnetization data is in good agreement with the 2D 
coupled dimer model. 



I. INTRODUCTION 



Quantum magnets with a dimerized singlet ground 
state are a fascinating research field with phenomena 
ranging from spin-Peierls transitions in antiferromagnetic 
spin 5 = 1/2 chains like CuGeOg 1 - or TiOCl 2 -^ to the 
Bose-Einstein condensation of magnons in external mag- 
netic fields as reported for the prototypical spin S = 1/2 
dimer systems BaCuSi2 0g^ and TlCuCl^ and also very 
recently for S^C^Os and BaaC^Og^r— . The singlet 
ground state and the excited triplet state of the dimcrs 
are separated by an energy gap due to the different ex- 
change energy between parallel and antiparallel aligned 
spins. This gap can be controlled by the external mag- 
netic field, which splits the triplet due to the Zeeman 
effect in such a way that the lower triplet level meets 
the singulet giving rise to the observed Bose-Einstein 
condensation* 1 - 1 -. This is a typical example of a quan- 
tum phase transition, i.e. it is driven by a non thermal 
control parameter, and, therefore, of high interest con- 
cerning the investigation of quantum critical behaviour. 

In the majority of dimer systems the dominating intra- 
dimer superexchange coupling is easily identified as the 
one mediated by the shared ligands of the magnetic 
ions in the structural dimer. In some systems, how- 
ever, this nearest-neighbour superexchange coupling is 
very weak and one has to consider more complex su- 
perexchange paths via several neighboring ions. It has 
been shown that the resulting exchange couplings can be 
of the same order of magnitude as the nearest-neighbour 
superexchange^ 2 ., with the spin-gapped CuTc20s - the 
compound under investigation in the present work - as 
one particular example 1 ^—. It is important to note 
that the influence of the nonmagnetic lone-pair tetrava- 
lent chalcogenide ion (Te 4+ ) can significantly change 
the superexchange interactions between the copper spins 
(Cu 2+ : S = 1/2) as seen for the case of related CuSe20.5 
(with Se 4+ ) which reportedly forms a spin-chain like 
magnetic structure with an antiferromagnetic transition 



at about 17 Kikil. 

In CuTe2 05, the situation turned out to be more 
complicated and the exchange geometry between the 
dimcrs (one- or two dimensional) remained unresolved 
so fa r 15 i 18 . Here we will show that our latest magnetiza- 
tion measurements are in favour of the two dimensional 
exchange geometry 



II. DESCRIPTION OF THE PROBLEM 

The crystalline structure of CuTe20s belongs to the 
space group P2i/c. The unit cell consists of four formula 
units, and its dimensions are a = 6.871 A, b = 9.322 
A, c = 7.602 A; the angle (3 between the a and c axes 
is 109.08°—. There are four nonequivalent positions of 
copper ions in the unit cell. Each of the copper ions is 
surrounded by six oxygen atoms forming a strongly dis- 
torted octahedron. Neighboring Cu pairs form structural 
CU2O10 dimer units, built from edge sharing octahedra 
which are rotated relative to each other and slightly mag- 
netically nonequivalent. The CU2O10 units are connected 
via Te 4+ ions (see Fig. [1]). Along the b axis the CU2O10 
units exhibit alternating rotation with respect to each 
other resulting in strongly magnetic nonequivalent Cu 
sites in neighboring units. 

The magnetic susceptibility of CuTe2 0s exhibits a 
broad maximum at about 57 K and an exponential de- 
crease to lower temperature o 14 i 1 5 . These data have been 
successfully described in terms of the model of a quasi- 
onc-dimensional antiferromagnetic spin chain with alter- 
nating exchange interactions, a modified Blcancy-Bowers 
model 1 ^ or a two-dimensional (2D) system of coupled 
dimers 1 ^. In all cases one obtains a leading exchange in- 
teraction of about 90 K, followed by a second exchange 
contribution of the same order of magnitude. The analy- 
sis of the susceptibility at the highest temperatures yields 
a Curie-Weiss temperature 9 = -41 K 15 i 19 . Electron 
spin resonance (ESR) studies evidenced the existence of 
nonequivalent magnetic Cu sites from the frequency de- 
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FIG. 1: Projection of the monoclinic lattice structure of 
CuTe205 with space group P2i/c on the fee-plane. The edge- 
sharing octahedra form CU2O10 units which are separated by 
Te ions. Ji, J4, and Ja indicate the most important exchange 
paths. 



evaluation of the data. 



III. EXPERIMENTAL DETAILS 

Single crystals of CuTe205 have been grown by a stan- 
dard chemical vapor phase method. Mixtures of high 
purity CuO (Alfa-Aesar, 99.995%) and Te0 2 (Acros, 
99.9995%) powders with a molar ratio higher than 1 : 
2 were sealed in quartz tubes with HBr as transport gas 
for the crystal growth. Then the ampoules were placed 
horizontally into a tubular two-zone furnace and heated 
very slowly by 50°C/h to 480°C for about 3 days un- 
der high vacuum. The optimum temperatures at the 
source and deposition zones for the growth of single crys- 
tals have been 580°C and 450°C, respectively. After two 
months, large blue-green plates with a maximum size of 
8 x 8 x 0.5 mm were obtained. X-ray powder diffraction of 
these crystals revealed the proper single-phase product. 

Magnetization measurements were performed in a DC- 
susceptometer (Oxford Instruments, Teslatron) in fields 
up to 14 T and in the temperature range 3-280 K. The 
temperature dependence of the specific heat of CUTC2O5 
was measured in a Physical Property Measurement Sys- 
tem PPMS (Quantum Design) between 1.8 K and 300 K. 



pendence of the linewidth dur to the anisotropic Zee- 
man effect and, therefore, suggested that strong superex- 
changc bonds should exist only between equivalent sites, 
while the supcrcxchangc between nonequivalent sites was 
estimated to be of the order of 0.5 K 15 i 20 i 21 . 

Starting from the crystal structure the isotropic ex- 
change interactions of nine different pairs of copper 
ions have been calculated by the Extended- Hiickcl-Tight- 
Binding (EHTB) method-^. With these calculations J 6 
was identified as the largest exchange interaction followed 
by J\ = 0.59 Jg resulting in alternating spin chains. From 
the fitting of the susceptibility data this estimate was ex- 
perimentally refined to J\ = 0.436J6, with J§ = 93.3 K. 
Later on the hopping integrals for the same pairs of 
copper ions in CuTe205 were obtained using density- 
functional calculations in combination with the Nth or- 
der muffin-tin orbital (NMTO) downfolding technique^ 
and LDA+U calculations^ 2 - and suggested that the sys- 
tem should be regarded as a two-dimensional coupled 
dimer system with dominant exchange coupling J4 fol- 
lowed by Je = O.28J4 and Ji = O.IU4. Again these es- 
timates were experimentally refined to Jg = 0.27 J4 and 
Ji = O.O7J4 with J4 = 94.2 K. Das and coworkers^ pro- 
posed to distinguish between the two magnetic models 
by their contribution to the magnetization and the spe- 
cific heat. In this work we compare both quantities the 
magnetization and the specific heat to these theoretical 
predictions. Note that both models postulate a strong 
exchange coupling between significantly nonequivalent 
copper sites, which seems to be in contradiction to the 
ESR results. We will comment on this problem after 



IV. RESULTS AND DISCUSSION 
A. Magnetization 

The predictions given by Das et alr^ for the magneti- 
zation of the 2D coupled dimer model were calculated for 
applied magnetic fields with strengths h/J = 0.2, 0.5, 1.0 
in units of the dominant exchange coupling J which cor- 
respond to absolute magnetic field values H = 12.7, 31.7, 
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FIG. 2: Temperature dependence of the magnetization mea- 
sured in 12.7 T of CuTe2 0s. The solid and dashed lines are 
calculations based on an alternating spin chain ( J§ = 93.3 K) 
and a 2D coupled dimer model ( J4 = 94.2 K) taken from Das 
et al. using the experimentally refined ratios of the exchange 
constants^. 
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63.4 T, respectively. Correspondingly, the magnetization 
for the alternating chain model has been calculated for 
H = 12.7 and 31.7 T for a direct comparison^. Quali- 
tatively both models yield a similar temperature depen- 
dence: starting from zero at T = due to the non mag- 
netic singlet ground state of the spin dimers, the magne- 
tization strongly increases with increasing temperature, 
when the dimers are breaking up, develops a maximum 
around the temperature corresponding to the antiferro- 
magnctic exchange coupling T max ~ J/fcrs, and finally 
decreases following a Curie- Weiss law. However, quan- 
titatively the absolute value of the maximum of the al- 
ternating spin-chain model (which is slightly shifted to 
higher temperatures) is about 25% below that of the 2D 
coupled dimer model. 

Only the lowest magnetic field is located in the range 
reachable in our standard experimental laboratory setup, 
while the higher fields have to be verified in special- 
ized high-field facilities. In Fig. [2] we show the data ob- 
tained at 12.7 T together with the predicted curves for 
the two models. The agreement of the experiment with 
the prediction for the 2D coupled dimer system is obvi- 
ously much better than with the one for the alternating 
spin chain based on ETHB calculations^. Therefore, the 
magnetization data clearly favour the 2D coupled dimer 
model as a realistic description of the magnetic structure 
of CuTe 2 5 . 



B. Specific Heat 

In Fig. Eta) we show the temperature dependence 
of the specific heat. No anomalies corresponding to 
phase transitions have been detected in agreement with 
reported susceptibility and ESR measurements^. As 
a result the specific heat increases monotonously with 
increasing temperature. Notably, at 300 K the spe- 
cific heat is still considerably lower than the expected 
high-temperature value 3i?s = 200 J/(mol K) for the 
phonon contribution given by the Dulong-Petit law indi- 
cating contributions to the phonon-density of states from 
higher-lying lattice modes. Here, R denotes the gas con- 
stant and s the number of atoms per formula unit. 

We assume that the total heat capacity originates from 
two different contributions, a lattice contribution C\ at t 
due to acoustic and optical phonons and a magnetic con- 
tribution C mag corresponding to the thermal population 
of excited dimer states. The magnetic contributions to 
the specific heat divided by temperature for a spin chain 
model C ait IT based on ETHB calculations^ and for a 2D 
coupled dimer system C2D/T, as predicted by Das and 
coworkers, are plotted in Fig. Etb) in comparison to the 
total specific heat C/T. Evidently, both magnetic con- 
tributions are small compared to the lattice contribution 
and a non-magnetic reference material is not available. 
Therefore, a straightforward method to unambiguously 
extract the magnetic contribution from the experimental 
data is difficult to realize. 
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FIG. 3: Temperature dependence of (a) the specific heat 
and (b) the specific heat divided by temperature (squares) 
together with theoretically suggested magnetic contributions 
for an alternating spin chain C a it/T (dashed) and a system of 
2D coupled dimers C 2 d/T (dotted) in CuTe 2 5 . The dash- 
dotted line in (a) indicates the high-temperature Dulong-Petit 
limit. The triangles in (b) denote the data after subtraction 
of the magnetic 2D coupled dimer contribution fitted by a 
pure lattice contribution (solid line). 



Thus, we decided to chose the reverse approach by sub- 
tracting the theoretically predicted magnetic contribu- 
tion for the 2D coupled dimer model, which is favoured 
by the magnetization measurements, from the total spe- 
cific heat and to analyze the residual lattice contribu- 
tion Ciatt- The resulting data were approximated fol- 
lowing standard procedures^ with a minimized set of 
fit parameters only using a sum of one isotropic Debyc 
(D) - accounting for the 3 acoustic phonon branches 
- and four isotropic Einstein terms (£1,2,3,4) - aver- 
aging the 3s — 3 = 21 optical phonon branches (fit- 
ting with less than four Einstein terms was not suffi- 
cient). For further reducing the number of free fit pa- 
rameters, the ratio between these terms was fixed to 
D : E x : E 2 : E 3 : E 4 = 1 : 1 : 2 : 2 : 2 to account 
for the 3s = 24 degrees of freedom per formula unit. For 
s = 8 atoms formula unit, the ratio between acoustical 
(Debye) and optical (Einstein) contributions is naturally 
fixed as 1:7. The weight distribution between the opti- 
cal contributions is chosen in such a way that the lowest 
Einstein mode is of equal weight with the Debye contri- 
bution, assuming one low-lying isotropic optical phonon 
branch (with 3 degrees of freedom), while the remaining 
degrees of freedom have been equally distributed between 
the higher Einstein modes. The resulting fit curve (solid 
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line in Fig. [3Jd) describes the data satisfactorily. For the 
respective Debye and Einstein temperatures we obtained 
Q D = 144 K, @ei = 151 K, Q E 2 = 308 K, 6^3 = 497 K, 
and Q E 4 = 1225 K. 

To check the fit results of the lattice contribution we 
compare the optical phonon frequencies detected by Ra- 
man spectroscopy in CuTe2 05 in the infrared region be- 
tween 100 < v < 900 cm" 1 to the Einstein tempera- 
tures obtained from the specific heat. The leading Ra- 
man peaks and corresponding characteristic tempera- 
tures (given in brackets) are found at v = 121 cm -1 
(174 K), 211 cm" 1 (304 K), 375 cm- 1 (540 K), 444 cm- 1 
(640 K), 745.5 cm" 1 (1074 K)2i. As one can see, the Ein- 
stein temperatures are in reasonable agreement with our 
fitting results, especially the low-lying optical mode at 
174 K corresponds to O^i = 151 K and also the ex- 
istence of high-frequency modes at 640 K and 1074 K 
agrees well with the fact that the Dulong-Petit value is 
approached only far above room temperature. 

V. DISCUSSION 

So far the present magnetization and specific-heat 
investigations favour the 2D coupled dimer model for 
CUTC2O5. The only open question remains concerning 
the frequency dependence of the ESR lincwidth which 
was explained in terms of an anisotropic Zccman effect 
requiring a very small exchange coupling between neigh- 
boring nonequivalent copper sites, whereas both theoreti- 
cal models derive a rather strong exchange between these 
sites. This dilemma can be probably resolved regarding 
very recent experimental findings. 

Current terahertz spectroscopic investigations of 
the singlet-triplet excitations in the dimerised low- 
temperature phase proof the importance of the 
Dzyaloshinsky-Moriya (DM) interaction in CuTe20g2£ 
which has not been taken into account in the previous 
analysis of the ESR lincwidth. Moreover, very recently 
the anisotropy as well as the temperature and frequency 
dependence of the ESR linewidth in the related uniform 
spin- 1/2- chain compound CuSe205 have been success- 
fully explained to result from the DM interaction^. Due 
to the theory of Oshikawa and Affleck, which models the 
spin relaxation in uniform spin-l/2-chains at low temper- 
ature, the lincwidth should increase with decreasing tem- 
perature and increasing frequency (i.e. resonance field) 
in the presence of a staggered field arising from the DM 
interactio n 27 ' 28 . It is especially important to note that 
the DM contribution to the linewidth increases propor- 
tional to the square of the applied frequency and cor- 
responding resonance field like the contribution of the 
anisotropic Zeeman effect. Thus, both contributions are 



difficult to separate. 

Looking at the cases of copper benzoate^, where Os- 
hikawa and Affleck first applied their theory, and also 
CuSe20^, it turns out that in both systems strong 
broadening with increasing frequency appears for the ex- 
ternal field along a certain crystal axis, while the broad- 
ening is significantly weaker for any perpendicular direc- 
tion. Despite the fact that CUT32O5 is not a uniform 
spin-l/2-chain, but a complex dimer system, the pro- 
nounced line broadening with increasing frequency for 
the magnetic field applied along the b axis compared to 
minor effects for the perpendicular orientations, suggests 
an analogous importance of the DM interaction for the 
spin relaxation in this compound. For a deeper anal- 
ysis additional theoretical effort is needed to check the 
applicability of the Oshikawa-Affleck theory to the case 
of CUTC2O5. At the moment we can state that the 
anisotropic Zccman effect alone is probably not enough to 
explain the frequency dependence of the linewidth and, 
hence, the exchange coupling between nonequivalent cop- 
per sites is not necessarily small. Insofar the seeming 
contradiction of the ESR data to the existing models is 
not mandatory. This finding supports our result that 
CuTe2C>5 is best described in terms of a 2D coupled dimer 
model. 



VI. CONCLUSION 

To summarize, the temperature dependent magneti- 
zation data of CuTe205 taken in an external magnetic 
field of 12.7 T agree well with the prediction for the 
2D-coupled dimer model and clearly deviate from the 
expectation for the alternating spin-chain model. Af- 
ter subtraction of the corresponding magnetic specific- 
heat contribution from the experimental heat-capacity 
data the residual specific heat can be satisfactorily fit- 
ted by phonon-contributions only, which are in line with 
the characteristic phonon energies obtained from Raman 
spectroscopy. Thus, the two-dimensional model is clearly 
favoured by the present experiments. 
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